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Abstract: The purpose of the paper is to construct a polyhedral Markov field in M 3 in 
analogy with the planar construction of the original Arak (1982) polygonal Markov field. 
We provide a dynamic construction of the process in terms of evolution of two-dimensional 
multi-edge systems tracing polyhedral boundaries of the field in three-dimensional time- 
space. We also give a general algorithm for simulating Gibbsian modifications of the con¬ 
structed polyhedral field. 


1 Introduction 


The notion of a consistent polygonal Markov field has first appeared in the seminal paper 
by Arak (1982) who constructed an isometry-invariant process with polygonal realisations 
in the plane, enjoying a two-dimensional germ Markov property [the conditional behaviour 
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of the field in an open bounded domain with piecewise smooth boundary depends on 
the exterior configuration only through arbitrarily close neighbourhoods of the boundary] 
and with one-dimensional sections coinciding in distribution with homogeneous Poisson 
point processes. An attractive feature of this process is that it admits several alternative 
equivalent representations. Two of them fall into the general Gibbsian framework and use, 
respectively, lines and points as the basic building blocks, see Arak & Surgailis (1989), 
Arak & Surgailis (1991) and Arak, Clifford & Surgailis (1993). The third representation 
of the Arak-Surgailis field is available in terms of equilibrium evolution of one-dimensional 
particle systems, tracing the polygonal realisations of the process in two-dimensional time- 
space, see the above papers. A dynamic representation in this spirit is in fact shared by 
a much richer class of planar polygonal Markov processes, as constructed and discussed in 
Arak & Surgailis (1989,1991) and Arak, Clifford & Surgailis (1993). 

The purpose of the present article is to construct an analogue of the Arak process in 
three-dimensional polyhedral setting. We provide a plane-based Gibbsian representation 
as well as a dynamic construction in terms of evolving multi-edge systems tracing polyhe¬ 
dral boundaries of the field in three-dimensional time-space Rxl 2 . The resulting field is 
Markovian, as directly follows by the Gibbsian representation. The construction is isom¬ 
etry invariant. We give also an explicit formula for the corresponding partition function. 
However, we are not able to explicitly characterise the equilibrium initial condition for the 
constructed dynamics for multi-edge systems and hence to explicitly construct a consistent 
version of our polyhedral Markov field. 

The results discussed above are presented Sections 2 and 3. They are complemented 
with Section 4 where we provide a Metropolis type simulation algorithm for a general class 
of Gibbsian modifications of the polyhedral Markov field constructed in this paper. Apart 
from its theoretical interest the purpose of developing this algorithm is related to the fact 
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that we anticipate possible applications of polyhedral Markov holds for volumetric image 
segmentation, in the spirit of our papers Kluszczynski et. ah (2004,2005), where similar 
algorithms in two-dimensional setting are used for planar image segmentation, see also 
Clifford & Middleton (1989) and Clifford & Nicholls (1994). 

2 Construction of the process 

2.1 Preliminaries 

In what follows we abuse the language granting the name polygonal face to any open 
connected subset of a plane in M 3 , with polygonal boundary. In particular, a face may well 
contain holes but it cannot split into several disconnected parts. For an open bounded 
and convex polyhedral domain D C M 3 we define the family of admissible polyhedral 
configurations in D by taking all collections 7 of planar polygonal faces contained in D 
such that 

(pi) the faces of 7 do not intersect, but they may share edges, 

(P2) each interior edge of 7 (contained in D) is shared by exactly two faces, 

(P3) each boundary edge of 7 (contained in dD) belongs to exactly one face, 

(P4) each internal vertex of 7 (contained in D) is shared by exactly three faces and has 
exactly three outgoing edges, 

(P5) each boundary vertex of 7 (contained in dD) is shared by exactly two faces and has 
exactly one outgoing internal edge and exactly two outgoing boundary edges, 

(P6) no two faces of 7 are coplanar. 

In other words, an admissible polyhedral configuration 7 corresponds to a family of disjoint 
closed polyhedral surfaces in D , possibly nested and chopped off by the boundary dD. We 
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write Ed{ 7 ) for the collection of edges of 7 contained in D (internal edges), Egu( 7 ) for the 
collection of edges of 7 contained in dD (boundary edges), and F D ( 7 ) for the collection of 
faces of 7 contained in D (only internal faces are considered). For a finite collection {vd ;}” =1 
of planes intersecting D we consider the family T £>({tUj}f =1 ) of all admissible polyhedral 
configurations 7 in D with the additional properties that 7 C (J ( l =1 w % and that 7 fl vd, is a 
single polygonal face of non-zero area, possibly with some zero-measure polygonal curves 
added, for each w i% i = 1,2, Let // be the usual isometry-invariant Haar-Lebesgue 
measure on the space V of all two-dimensional planes in R 3 . One possible construction of 
H goes by identifying a plane vd with the pair (u, p) 6 §2 x R + , with § 2 standing for the 
unit sphere in R 3 and where u is the direction of the vector orthogonal to vd and joining 
it to the origin, while p is the distance between vd and the origin. In these terms p arises 
by endowing § 2 x R + with the product of the usual surface measure on the sphere and the 
Lebesgue measure on R + . We write II for the Poisson plane process with intensity measure 
/i and for each bounded domain D we let IIstand for the restriction of II to the family 
Vd of planes intersecting D. The following properties of the plane process II will be of 
use in the sequel. Here and throughout, to avoid possible confusion we use d 2 x and d 3 x 
rather than dx to denote integration respectively w.r.t. the 2- and 3-dimensional Lebesgue 
measure, thus explicitly indicating the dimensionality of the integration variable. 

Proposition 1 With the above notation we have: 

(11) The intersection of n with a given straight line l in R 3 is a homogeneous Poisson 
point process of intensity 7 r, 

(12) The intersection of n with a given plane vd is a Poisson process of lines l = 1(f), r) 
in vd with intensity measure | df>dr , where r stands for the distance between l and a 
certain fixed point 0 in vd, while 0 is the angle between some fixed line in vd and the 
vector in vd orthogonal to l and joining it to 0, 
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(13) For a given plane w and x G w the probability that two planes of II meet at [x, x + 

3 

d 2 x] C w is ^d 2 x. 

(14) For x G M 3 the probability that three planes of II meet at [x, x + d 5 x] C M 3 is ^d 3 x. 

Indeed, the assertions (II) and (12) follow immediately by (3.29T) in Miles (1971) because 
our process II coincides with B(2n, 2, 3) there. Further, (14) is a direct consequence ol ( 6 - 
1-9’) in Matheron (1975) with a = 7 t /2 there. Finally, (13) follows from (12) and (6-1-9’) 
ibidem with a — 7t/ 2 there. 

2.2 Gibbsian representation 


For a given admissible polyhedral configuration 7 we consider its energy 


® d ( i ) = ~ l Z ( e )IK( e ) + Tf Area (/) + y Vol(.D), ( 1 ) 

eS-E D ( 7 ) f&F D {l) 

where Z(e) is the solid convex angle between the planes determined by the faces meeting 
at an internal edge e, £(e) is the length of the edge e, Area(/) stands for the area of 
the polygonal face / while Vol(Z7) denotes the volume of the domain D. We define our 
three-dimensional polyhedral process Ad in D as the Gibbsian modification of the process 
induced on by IId with the Hamiltonian 7 1 —> 4 * 1 )( 7 ). To be more specific, we put 

E E 7e r D (n D )n^ ex P(-M7)) 


P (Ad EF) = 


( 2 ) 


E E 7G r D( n D ) ex P(-3fo(7)) 
for all T C T d Borel measurable, say, with respect to the standard Hausdorff distance 
topology. The finiteness of the partition function E X) 7 er D (ii D ) ex:p(—above will be 
established in Section 2.4, Corollary 1. 


2.3 Dynamic representation 

The availability of a dynamic description in terms of equilibrium time-space evolution of 
one-dimensional particle systems is a crucial feature of the Arak-Surgailis process and it 
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underlies most of the techniques successfully used to derive explicit expressions for var¬ 
ious characteristics of that model. It is therefore of particular importance to show that 
our three-dimensional model Ad admits an analogous dynamic representation in terms 
of evolution of two-dimensional edge systems tracing the polyhedral boundaries of Ad hr 
three-dimensional time-space. 

Note that a planar section of a single polygonal face as defined above coincides, for 
//-almost all planes in V, with a finite collection of disjoint colinear closed segments. Moti¬ 
vated by this observation we shall use the name of multi-edge for each such collection. For 
the purposes of this section we shall represent M 3 as three-dimensional time space MxR 2 
with two-dimensional spatial component; the time coordinate will be consistently denoted 
by t. All the multi-edges considered below will arise as sections of polygonal faces in ixl 2 
with purely spatial planes Wt {£} x M 2 . In this setting, a non-empty multi-edge given 
as the intersection £ := f D w t for some polygonal face /, / (Z w tj can be assigned in a 
natural way its two-dimensional velocity vector lying in zu t and defined to be the spatial 
component (projection onto the spatial plane zu t ) of the unit normal to £ lying in the 
plane pl[/] of / and contained in the future half-space w t + [t, +oo) x M 2 , divided by 
the length of the spatial component of the unit normal to pl[/]. With g standing for the 
unique plane pl[/] containing the face / we denote by v[(f\ the velocity vector constructed 
above (note that the construction does not depend on t). It is now easily verified that 
when observing the time evolution of the spatial plane w t , we see the polygonal face / 
mapped onto the time evolution of the multi-edge wt D / moving in time with the velocity 
vector u[pl(/)]. Consequently, each finite family of disjoint and possibly nested bounded 
polyhedra is also mapped into the time evolution of the corresponding finite multi-edge 
systems. In course of the time evolution of multi-edges their components evolve as well, 
continuously updating their lengths according to the rule that no two edge segments can 
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intersect each other except when sharing a vertex, whence the endpoints of a given edge 
segment are determined as the meeting points with other multi-edges. This rule may lead 
to extinction of certain edge segments as their length reaches 0 but, on the other hand, 
certain segments may also get split into disjoint parts. A multi-edge dies if all its edges go 
extinct. 

As mentioned above, we shall consider our polyhedral random field defined in a certain 
open bounded and convex polyhedral domain D C I x l 2 . The presence of the boundary 
has its effect on the dynamics, which has to be taken into account in our construction. 
Namely, at each time t El we observe a collection of boundary multi-edges in w t r)dD and 
for each internal multi-edge in Dr\zu t the meeting point with the boundary marks the end of 
its appropriate segment. The boundary multi-edges can be assigned their velocity vectors 
on equal rights with the internal edges, ft should be noted that a boundary multi-edge has 
always exactly one segment due to the convexity assumption imposed on D. 

Apart from the multi-edge extinction, the evolution of a multi-edge system, correspond¬ 
ing to a finite collection of disjoint and possibly nested polyhedral surfaces, comprises also 
multi-edge birth events. These can be divided into several groups 

(IT) Infinitesimal triangle birth in a time-space point x G D n w t which does not lie on 
any internal or boundary multi-edge in w t . At the point x three pairwise non-parallel 
multi-edges e±, e 2 , e% are born, moving with three different velocities Vi, V 2 , # 3 . 

(IA) Infinitesimal angle birth in a time-space point x lying on an internal or boundary 
edge ei moving with a certain velocity v\. At the point x two multi-edges e 2 ,e 3 are 
born, moving with velocities fq, iq and such that {ei, e 2 , e^} are pairwise non-parallel. 

(IE) Infinitesimal edge birth at a vertex in a time-space point x lying at the intersection 
of two different (either internal or boundary) edges ei,e 2 C (T> U iID) n w t , moving 
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with velocities V\ and w 2 respectively. At the point x a new multi-edge 63 is born, 
moving with velocity v 3 and such that {ei,e 2 ,e 3 } are pairwise non-parallel. 

Observe that for certain choices of new-born multi-edges and their velocities we end up with 
unstable configurations where new edges die immediately upon their birth. This happens 
e.g. for an infinitesimal triangle with all its edges moving inward. We say that a multi-edge 
e is stably born if it survives a positive amount of time rather than going extinct just upon 
its birth. The stability conditions are easily determined for all groups of birth events. 

(Stability for IT) Denote by n; the unit normal to e* (in w t ) pointing outward the 
infinitesimal triangle and let v [vj , v k ] be the velocity vector (in wt) for the intersection 
point of ej and e k , easily checked to be avj + bv k , where 


a 


(O'Dh) 

{Vj,V k ) 

-1 

f lof 1 

b 


. (Vk,Vj) 

( Vk,V k ) _ 


L i«*i 2 J 


The birth is stable iff (v),rij) > (v[uj, v k \, n.j) for all i j k i. It can be easily 
shown that it is enough to verify this condition just for one fixed choice of i,j, k. 

(Stability for IA) Let ni be the unit normal to e\ in Wt pointing to the side of e± 
opposite to where the new angle is born and write n 2 ,n 3 for the unit normals to 
e 2 , e 3 respectively, pointing outward the new-born convex angle. Once again, the 
birth is stable iff (v t , n, : ) > (v[vj, v k ], n*) for all 1 j k i. 

(Stability for IE) Let n 3 be the unit normal to e 3 pointing inward the convex angle 
between e\ and e 2 . Then the birth is stable iff (n 3 ,n 3 ) > (v[vi, n 2 ], n 3 ). 

In geometrical terms the stability condition for an mfinitesimal triangle ei,e 2 ,e 3 means 
that e, is able to escape in the direction pointed by n, from the intersection point of ej and 
e^, whose velocity component in the direction n* is smaller than (v t , n*). Clearly, otherwise 



would be destroyed by other edges immediately upon its birth. The remaining stability 
conditions admit their geometrical interpretations along the same lines. 

In the sequel, we will refer to all stable boundary birth events (involving at least one 
boundary edge) of either (IA) or (IE) type as to entry events. Note that (IA) or (IE) 
events involving internal edges only are not considered to be entry events and neither are 
unstable birth events. For a given admissible polyhedral configuration 7 € T n we write 
Entry (7) for the collection of entry events it determines. 

Below, we construct a random dynamics on multi-edge systems conditionally on the 
collection £ of all entry events assumed to be given. This stands in contrast to the two- 
dimensional Arak-Surgailis construction, where boundary birth events were also governed 
by explicit random dynamics, and this is due to the fact that we are unable to identify the 
explicit equilibrium distribution for entry events in the three-dimensional case. 

The random dynamics of internal birth events is given by the following rules, with t 
standing for the time coordinate 

(Dynamics for IT) The infinitesimal triangle birth time-space sites are chosen according 

4 

to a homogeneous Poisson point process of intensity in D. The directions and 
velocities of the new-born edges are chosen according to the distribution of a typical 
stable vertex angle at a meeting point of three planes in II, where the adjective stable 
means here that out of the 8 vertex angles at the triple intersection point we choose 
the unique one which gives rise to a stable infinitesimal triangle birth in the sense of 
(Stability for IT), 

(Dynamics for IA) The infinitesimal angle birth time-space sites on a face / are chosen 
with intensity ^a /1 + \v[pl(f)}\ 2 d£dt with d£ standing for the edge length element 
on / fl w t (note that this intensity coincides with times the area element on 
/). The directions and velocities of the new-born edges are chosen according to the 
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distribution of a typical stable vertex angle at a meeting point of pl(/) with two other 
planes of II, where the adjective stable means here that out of the 8 vertex angles 
at the triple intersection point we choose the unique one which gives rise to a stable 
infinitesimal angle birth in the sense of (Stability for IA), 

(Dynamics for IE) An infinitesimal edge birth site at a vertex x e try, at the intersection 
of two edges e\ and e 2 lying respectively on faces f\ and f 2 and moving with respective 
velocities v\ and v 2 , arises with intensity 

7Ta/1 + | v[vi,v 2 \\ 2 d 2 xdt. (3) 

Note that this intensity coincides with n times the length element on the line traced by 
the time-space trajectory of the intersection point e\ fle 2 . The direction and velocity 
of the new-born edge are determined according to the distribution of a typical plane 
w of II intersecting the line traced in time-space by e± De 2 . Unlike in previous cases, 
the stability of this birth event in the sense of (Stability for IE) is not guaranteed 
here and in fact it is easily seen that it holds iff the normal to w pointing in time- 
space to the future is not contained in the solid angle between the faces f\ and f 2 
[i.e. the solid angle between the half-planes of pl[/i] and pl[/ 2 ] containing e\ and e 2 
respectively and meeting along pl[/i] fl pl[/ 2 ]], denote this angle by Z(/i, f 2 ) for use 
below. The non-stable birth events resulting from the above rule have no effect on 
the dynamics. Thus, alternatively we can produce the birth events with intensity 

, i r -fP7(2vr-|Z(/i,/ 2 )|) 

7T \J 1 + |v[ui,u 2 J| 2 --- d 2 xdt (4) 

Z7T 

and let the direction and velocity of the new-born edge be determined according to 
the distribution of a typical plane w of II intersecting the line traced in time-space by 
eifle 2 , conditioned on yielding a stable infinitesimal edge birth event. In analogy with 
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the interpretation of (3) above, the intensity (4) coincides with |(27 t — |Z(/i,/ 2 )|) 
times the length element on the line pl[/i] npl[/ 2 ] traced by the time-space trajectory 
of the intersection point e\ fl e 2 . 

To make the above construction fully explicit, we note that the distribution of a typical 
vertex angle between three planes of fl is oc |(ni,n 2 x n 3 ) |da(ni)da(n 2 )<icr(n3), where 
rij G § 2 , ° is the surface measure on § 2 and the planes i = 1,2,3 creating the an¬ 
gle are respectively chosen orthogonal to rq, i = 1,2,3; indeed, this follows by a minor 
modification of Theorem 3 in Calka (2001), specialised for d — 3. Recall that the scalar 
product |(n 1 ,n 2 x n 3 )| coincides with the volume of the parallelepiped spanned by the 
vectors ni,n 2 ,ri 3 . Observe that the knowledge of the typical vertex angle distribution for 
three planes provides full knowledge of the corresponding laws with one or two planes fixed, 
as respectively required for (IA) and (IE). 

The polyhedral process in D resulting from the above dynamic construction in presence 
of a collection £ of entry events will be denoted in the sequel by Ad\s- 

2.4 Equivalence of representations 

We now proceed to showing that both the Gibbsian representation (2) and the dynamic 
construction of Subsection 2.3 yield, in a sense to be specified below, the same polyhedral 
field. To this end, for a given collection £ of entry events in D we put V D \£ := {7 G 
T o | Entry( 7 ) = £}. Further, for a collection {zu j }” =1 of pairwise non-parallel planes in M 3 
we write Td\s ((tu}" =1 ) to denote the family of all admissible polyhedral configurations in 
V D \£ with the additional property that (pl(/) | / G F D ( 7) \ F[£]} coincides with (w}" =1 , 
where F[S] is the collection of planes arising in entry events from £. Define k(D) to be the 
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measure /i of the set of planes intersecting D. Note that, D being convex, we have 

K i D ) ■= ufa,™ n D ± 0} = ^ Y l Z ( e )l^( e )> ( 5 ) 

e£E{D) 

see (4.2.30), (4.5.9) and (4.5.10) in Schneider (1993), whence n(D) is proportional to the 
generalised integral mean curvature functional of D, see ibidem. 

The crucial observation in this Subsection is that, in view of Proposition 1, the form of 
the dynamic rules (Dynamics for IT,IA,IE) implies that, given a collection £ of entry 
events in D, we have for all 7 e T £>\£ 

P (A D \e e d'y) = exp(K(D))l {7e r D|£ (n D )} exp(-<f> D (7))d£(n D ) (6) 

with jC(X) standing for the law of a random object X. Indeed, for 7 e £d\£, the proba¬ 
bility differential [dC{AD\e)\{l) factorises into the product of the probability that all faces 
/i, - - - Jk of Ffj ( 7 ) \ F[£\ were born in course of the evolution of the multi-edge system, 
which is seen to be E[/=i [ 1 P i[/dn^^ 0 ^(pl[/i] )] = exp(K(D))[d£(A D )]({pl[/ 1 ],... ,pl [/*]}) by 
comparing the dynamic rules (Dynamics for IT,I A,IE) with Proposition 1, times the 
probability that no other faces were born in the evolution giving rise to 7, which is in its 
turn evaluated to exp(—$£>( 7 )) in view of the form of the dynamic rules. This last claim 
is easily verified by noting that 

• by (Dynamics for IT), exp ^ Vol(D)^ is the probability that no extra faces 
were born in (IT) birth events, 

• by (Dynamics for IA), exp (^ — x YlfeF D ( 7 ) Area(/) j is the probability that no extra 
faces were born in (IA) birth events, 

• by (Dynamics for IE), exp | X^ e e£' D ( 7 )( 2 7r — |A(e)|)^(e)j is the probability that 
no extra faces were born in (IE) birth events, see the discussion following (4). 
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As an immediate consequence of ( 6 ) we come to 


Theorem 1 Given a collection £ of entry events in D, the random polyhedral field Ad\s 
coincides in distribution with the Gibbs polyhedral field A* D ^ £ given by 

E E 7e r D|f( n D )n^ ex P(-Mx)) 


P (A% e e F) = 


E E 7G r 0|f( n D ) ex P(-^(7)) 


for each T C T £>\s Borel measurable with respect to the usual Hausdorff topology. Moreover, 
for each £ we have 


E Y exp(—$£>( 7 )) = exp(-/c(D)). (7) 

7er D | £ (n D ) 

It should be noted at this point that the particular form of the expression (7) for the 
partition function might seem to stand in an unexpected contrast to the two-dimensional 
formula (4.6) in Arak & Surgailis (1989). However, this difference is a matter of choice of 
the reference measure, which is in our case the normalised law of the Poisson plane process, 
while Arak & Surgailis (1989) use instead.. Moreover, we introduce the volume 

order term ^ Vol(-D) to our energy function, which has no equivalent in that paper. We 
find these choices preferable for the presentation of our setting, as leading to simpler 
formulae. However, should we use the reference measure and Hamiltonian analogous to 
those of that paper, our partition function in (7) would evaluate to exp Vol(T>)^ . The 
absence of k(D) and a surface area order term in the exponent in this partition function 
expression is due to the fact that we condition on fixed collection of entry events rather 
than randomising it as done in Arak & Surgailis (1989). 

As an obvious conclusion from Theorem 1 we get that the original polyhedral held Ad 
as defined by ( 2 ) admits the dynamic representation with the collection of entry events 
distributed according to Entry (Ad). 
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Corollary 1 We have 

C{A d ) = j C(A D \e)d[C(EntTy(A D ))]{£) 

and 

£{Ad\e) — £{Ad \ Entry(A-D) = £). 

Moreover, 

E exp(-$ D ( 7 )) = exp (-k(D)). 

7er D (n D ) 

The use of this corollary is limited by the fact that we do not know the distribution of 
Entry (,A d ). 

3 Properties of the process 

We argue first that the polyhedral fields Ad and Ad\s exhibit a 3-dimensional germ-Markov 
property in the sense specified in Corollary 2 below. For a smooth closed simple (non-self- 
intersecting) surface a in a bounded and convex polyhedral domain D, by the trace of a 
polyhedral configuration 7 on a, denoted in the sequel by 7 A a, we mean the intersection 
7 flcr together with the directions of normals to the face planes at intersection. This concept 
can be formalised in various compatible ways, yet we keep the above informal definition 
in hope that it does not lead to any ambiguities while allowing us to avoid unnecessary 
technicalities. For convenience we assume that no face of 7 is tangent to c, which can 
be ensured with probability 1 in view of the smoothness of a. In view of the Gibbsian 
representation (2) and by Theorem 1 we easily conclude that 

Corollary 2 For each a as above there exists a stochastic kernel Ai n ta(-\^) such that, with 
tt standing for a trace on a, 

Ant a (Ad\Ad A cr = l9) — Ant a {Ad\e\Ad\£ A <J = $) = Ant a ( ' | A (8) 
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for all bounded open and convex polyhedral domains D D Int a and for each collection £ of 
entry events in D, where C\ nt(7 {X) denotes the law of a random field X restricted to Inter 
(the interior of a). 


Indeed, it is easily seen that we have for measurable G C T Into .p 

4 „„ E E 76 r lat . ll , ( n In ,„ l nGe x P(- , l > iM«(r)) 

( 9 ) 

EE 7 e r IntCT|l ,(n IntCT )exp(-$i ntff ( 7 )) 

where T Int(T p := {7 G Tinto- | 7 Aer = d} and, for a collection {wj }" =1 of planes in M 3 hitting 
Int cr, r Int( 7 p({wj}” =1 ) denotes the family of all polyhedral configurations 7 G T Int(T p for 
which the set {pi[/] | / G Fp( 7 ), / fl a = 0} coincides with {tZ7 ,}” =1 • 

To proceed, consider the family T R 3 of whole-space admissible polyhedral configurations, 
determined by (PI), (P2), (P4) and (P6) ((P3) and (P5) are meaningless in this 
context) and by the requirement of local finiteness (any bounded set is hit by at most a 
finite number of faces). It is natural to define the family G(A) of infinite volume Gibbs 
measures (thermodynamic limits) for A as the collection of all probability measures on T R 3 
with the accordingly distributed random element A satisfying 


Ant a (-4-1-4. A a = -d) = 4int<r(-|0) (10) 

for a ranging through the collection of all bounded smooth simple closed surfaces in M 3 . 
In addition, we shall consider the family Q T {A) of isometry invariant measures in Q(A). 
We believe that, in analogy with the results in Section 3 in Schreiber (2005), it should be 
possible to show that G T (A) ^ 0 by using an appropriate relative compactness argument. 
Moreover, we conjecture that the uniqueness of the isometry invariant thermodynamic 
limit Mk 3 , as well as its coincidence with the polyhedral process traced by infinite-volume 
equlibrium evolution of the multi-edge system as discussed above, could possibly be es¬ 
tablished following the lines of Schreiber (2004a), where this is done for two-dimensional 
polygonal fields admitting dynamic representation. On the other hand, in analogy with 
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the two-dimensional setting, we do not expect that 4 R 3 be the unique element of 6(A), 
see the discussion closing Section 3 in Schreiber (2005). We are working on this conjecture 
at present, yet we are unable to provide their formal proofs at the current stage of our 
research. However, should these conjectures hold true as stated, initiating the multi-edge 
system dynamics in a domain D with the collection of entry events Entry(„4 R 3 D.D), for 4 R 3 
denoting the unique thermodynamic limit, would result in a consistent family of polyhedral 
fields, as constructed by Arak & Surgailis (1989) for the polygonal setting. A further es¬ 
sential task would be to provide a feasible description of the entry process Entry (4*3 fl D), 
which we anticipate to be of a rather complicated nature. 

4 Birth site birth and death dynamics for simulating 
polyhedral fields 

The purpose of the current section is to construct, much along the lines of Schreiber (2005), 
Section 2.1, a random dynamics on the space T D of admissible polyhedral configurations 
which leaves invariant the law of Ad\s , where the collection £ of entry events is to remain 
fixed throughout the section. This will allow us later to provide modifications of this dy¬ 
namics suitable for simulation of Gibbsian modifications of Ad\s- The purpose of developing 
this algorithm is its envisioned application, as a component of suitable simulated anneal¬ 
ing techniques, to volumetric image segmentation along the lines of our previous papers 
Kluszczynski et al. (2004,2005) where we considered the corresponding two-dimensional 
problem. 

In the sequel, particular care is needed to distinguish between the notion of time con¬ 
sidered in the dynamic representation of the polygonal held Ad\s given in Subsection 2.3 
above, and the notion of time to be introduced for the random dynamics on T o\£ con¬ 
structed below. To make this distinction clear we shall refer to the former as to the 


16 



representation time (r-time for short) and shall keep for it the notation t, while the latter 
will be called the simulation time (s-time for short) and will be consequently denoted by s 
in the sequel. 

It is convenient for onr exposition below to perceive each individual infinitesimal triangle 
birth site ((IT)-birth site) in the dynamic representation, see (IT), as coming with an 
associated random number generator, represented for instance as an infinite sequence of 
i.i.d. random variables uniformly distributed on [ 0 , 1 ] and used to determine the subsequent 
moments and angles/velocities for critical events (IA) and (IE) involving multi-edges 
resulting from the considered (IT) birth event. In other words, each (IT)-birth site is 
assumed to carry a package enclosing all randomness the resulting multi-edges may possibly 
encounter during their evolution, and the above is just one technical possibility of how this 
can be achieved. We shall use the name of a birth package for an infinitesimal triangle birth 
site with such a random number generator attached. In these terms, it is now easily seen 
that the polyhedral configuration obtained in course of the dynamic construction depends 
deterministically on the underlying collection of birth packages. 

Consider a polyhedral configuration 7 6 Yo\e and a new infinitesimal triangle birth site 
x 0 G D not yet present in 7, extended to a birth package in the standard way as discussed 
above. Adding this birth package to the collection of birth packages determining 7 and 
keeping the evolution rules of the dynamic representation (Dynamics for IT,IA,IE) 
results in a new configuration to be denoted by 7 © x 0 - Likewise, removing an (IT) -birth 
site X\ from a configuration 7 in which it was present yields a new polyhedral configuration 
7 © Xi- 

Taking into account that the collection of the (IT)-birth sites for Ao\e is chosen ac¬ 
cording to the Poisson point process with intensity ^ as specified in (Dynamics for IT), 
we easily see that the law of Ad\s is invariant with respect to the following pure-jump 
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Markovian (IT)-birth site birth and death dynamics on rz?|£r, further denoted by (BS), 
with standing for the state at time s and with y 0 ^ Fd\£ 

(BS:birth) With intensity ^ds set 7 s+ a s 7 s © x, 

(BS:death) For each (IT) -birth site in % with intensity 1 • ds set 7 s+ d s 7s © x, 

otherwise keep 7 s+ d s = %■ 

In fact, more can be stated, see also Proposition 1 in Schreiber (2005) 

Theorem 2 The distribution of the polygonal field Ad\s is the unique invariant law of 
the dynamics given by (BS:birth), and (BS:death). The resulting stationary process is 
reversible. Moreover, for any initial distribution of 70 concentrated on T D \ £ the laws of the 
random polygonal fields % converge in total variation to the law of Ad\s as s — > 00 . 

While the invariance was discussed above and the reversibility is clear, the uniqueness 
and convergence statements in the above theorem require a short justification. They both 
follow by the observation that, in finite volume, regardless of the initial state, the process 
7 S spends a non-null fraction of time in the state where no polyhedral faces other than 
those arising in £ are present. Indeed, this observation allows us to conclude the required 
uniqueness and convergence by a standard coupling argument. 

Dynamics for Gibbsian modifications Assume that a Hamiltonian (energy function) 
Tt is defined on the space Td of admissible polyhedral configurations and it satisfies 

H{i) > ~A\o\(D) -B (11) 

for some positive constants A, B. Then it is clear that the partition function 

Z D [h\ := Eexp(—7-f(M D |£-)) 
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is finite. Consequently, the corresponding Gibbsian modification A^ £ can be considered 
with 


d£(A™\ £ ) 

(1C(Ad \£) 


exp(—77(y)) 
Z d [H] 


7 e r D \ e . 


Consider the following modification of the basic (BS) dynamics: 


( 12 ) 


(BS[H]:birth) With intensity ^ ds propose the update <5 := 7 s ©x. Then, with probability 
min(1, exp(77(y s ) — 77(5))) accept this update, putting 7 s+ d s ■= 5, otherwise keep 
Ts+ds • 


(BS[H]:death) For each (IT) -birth site in % with intensity 1 • ds set 5 := 7 s Qx. Then, 
with probability min(l, exp( 77 ( 7 s ) — 77(5))) put 7 s+ d s '■= 8, otherwise keep 7 s+ d s '■= %, 


In other words, the original dynamics (BS) is used in the standard way to propose a new 
configuration 8, which is then accepted with probability min(l, exp( 77 ( 7 s ) —77(5))) and 
rejected otherwise. As a direct consequence of Theorem 2 and by a standard check of the 
detailed balance condition we get 

Theorem 3 The distribution of the polyhedral field A^ e is the unique invariant law of the 
dynamics given by (BS[H]:birth) and (BS[H]:death). The resulting stationary process 
is reversible. Moreover, for any initial distribution of 70 concentrated on T d\ e the laws of 
the random polyhedral fields 7 S converge in total variation to the law of A™> £ as s —> 00 . 
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